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Abstract 

A general approach for construction of quark and lepton mass matrices is formulated. 
The hierarchy of quarks and charged leptons ("electrons") is large, it leads using the 
experimental values of mixing angles to the hierarchical mass matrix slightly deviating 
from one's suggested earlier by Stech and including naturally the CP-phase. 

The same method based on the rotation of generation numbers in the diagonal mass 
matrix is used in the electron-neutrino sector of theory, where neutrino mass matrix is 
determined by the Majorano see-saw approach. The hierarchy of neutrino masses, much 
smaller than for quarks, was used including all existing (even preliminary) experimental 
data on neutrinos mixing. 

The leptonic mass matrix found in this way includes not known value of the leptonic 
CP-phase. It leads to a large v^v T oscillations and suppresses the v e v T and also v e Vn, 
oscillations. The explicit expressions for the probabilities of neutrino oscillation were 
obtained in order to specify the role of leptonic CP-phase. The value of time reversal 
effect (proportional to sin 5') was found to be small ~ 1%. However, a dependence of the 
values of ^f^fei/,- transition probabilities, averaged over oscillations, on the leptonic 
CP-phase has found to be not small - of order of tens percent. 
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1. Introduction 



Serious efforts have been invested recently in the natural understanding of experimental 
results on neutrino oscillations. They have shown that neutrinos of three generations have, 
perhaps, non vanishing small masses. The heaviest of them, the neutrino of the third genera- 
tion, seems to have a mass of the order of ~ (l/20)eV and, as the Super Kamiokande data on 
atmospheric neutrinos show, has the maximal possible mixing with the neutrino of the second 
generation, and may be, also not too small mixing with that of the first generation. This was 
not expected a priori, since all similar mixing angles of quarks are small. 

It is the challenge of modern particle physics to include naturally these results into the 
framework of Grand Unification Theory together with the data on quark masses and mixing 
angles. For the quarks these angles are small and are known already [1]. 

We begin this paper by reminding the well-known picture of masses, CP-phase and mixings 
for the quark sector of a theory. A general method will be developed which allows one to 
construct consistently the 3x3 mass-matrix and the CKM mixing matrix for quarks. The 
same general approach will be used later for the electron-neutrino sector of a theory. 

Let us consider, as a useful introduction to a consistent theory of quark and lepton masses 
and mixings, a simple phenomenological approach suggested by B. Stech [2]. He has noticed 
the following quark and charged lepton ("electrons") mass hierarchies: 



m t : m c : m u ~ 1 : o : a 



m b : m s : m d ~ 1 : —a : 8a" 



m T : : m e ~ 1 : a : -a 3 , (1) 



with a very small a 2 ~ 1/300, a ~ 0.058. He has also introduced the following mass matrices 
which reproduced approximately the masses of all the quarks as well as theirs mixing angles: 
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Here 77 = e lS represents the CP violating phase 5 in the quark sector, while the values of 
the constants a d ~ 2, a e ~ correspond to the best fit of all masses of quarks and electrons 
(their central values, see below). 

The diagonalization of the matrices (2) and (3) by means of an unitary matrix U a : 



Mdiag = (j a M a (j+ 



a = u,d, e, 



(4) 



1 



(where tfjj^ = UaU a = 1), reproduces simultaneously both the experimental (central) val- 
ues of running masses of quarks and electrons and all quark mixing angles (their sines): 
S12 = sin i?f 2 > S23 = sin ^23) s i3 — si n $13 m the CKM matrix: 



VcKM — U u Ud 



( C12C13 S12C13 si 3 e l5 ' \ 

-S12C13 - Ci 2 S23Sl 3 e i ' 5 ' C12C23 - Si2S2 3 Si 3 e i ' 5 ' S 23 Ci 3 
V S12S23 — C 12 C23Si3e" S — C12S23 — Si2C2 3 Si3e*' 5 C23C13 j 



(5) 



Here, according to experimental data [1,6], one has 0: 



S12 ~ sin ^ 2 ~ 0.221 ± 0.004, s 23 = sin 1^3 = 0.039 ± 0.003, 
S13 = sin^? 3 ~ si 2 s 23 |i?| = 0.0032 ± 0.0014, \R\ = s 13 /s 12 s 23 ~ 0.38 ± 0.19 (6) 
and 5 = ix/2 ± 7r/4. 



The factors l/v 7 ^, 1/2, 2, ^3/2, etc. in Eq.(2),(3) were adjusted so as to reproduce all these 
mixing angles and all the masses of quarks and leptons. As has been mentioned, for quarks 
the mixing angles are small, i.e. s 13 s 2 3 <C si 2 <^ 1 and therefore the cosines C13, c 2 3, ci 2 of 
them can be put approximately equal to unity. 

The major ideas of Stech's approach were: 1) to consider the running masses of all quarks 
at the same scale (instead of the pole ones) and 2) to obtain these masses and mixing angles 
in terms of few parameters. However, the matrices (2) and (3) have been written phenomeno- 
logically, just by hand. 

Below we develop an approach to quark and lepton masses and mixings through the 
following steps: 

a) A general method is suggested which allows one to construct the quark mass matrices 
consistently (without using Eq.(2)) in terms of their masses and mixing angles (6). The 
resulting mass matrices coincide with those given by Eq.(2). They are diagonalized in an 
analytical form of decomposition in powers of small parameter a and the quark CKM matrix 
is reproduced analytically. 

b) The same approach is applied to the reconstruction the neutrino mass matrices M v . 
This is performed using the Eq.(4) and the unitary matrix U v . Preliminary results on masses 
of neutrinos and their mixing angles, extracted from neutrino oscillation experiments have 
been used. 

c) The leptonic CKM matrix Vq KM = U^U e containing the leptonic CP-phase 5' is in- 
troduced. It leads to z/ M z/ r oscillation of the type observed at Super-Kamiokande and also to 
some suppression of v^v^ oscillation. 



lr The upper index q is used to emphasize that the angles i?y = i??-, or = sini!^- determine just by the 
quarks and not by the leptonic mixings. The corresponding leptonic angles and their sines are denoted below 
by i^ - and by = sini?^ without any upper index (see sect. 3 and 4). The leptonic CP-phase is denoted 
(in sect.3,4) by 5'. 
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d) The expressions for the three type neutrino oscillations are specified and some realistic 
numerical examples of oscillation dependence on the leptonic CP-phase are presented. 

Finishing this section let us note that the electron and quark masses emerge from the 
usual Yukawa interaction: 

L Yuk = ViUih^qj + (pxditijqj + ^xe^lj , 

where our mass matrices (2) and (3) are: (M u ) = {f2jh u , (M d ) = {tpi)h d , (M e ) = (ipi)h e 
and ((fix) = v cos f3, (</? 2 ) = v sin/3 are V.E.V of SUSY two neutral CP-even Higgs fields 
(v = 174GeV, and tg (3 are two well known SUSY parameters) and h l J ~ h d ~ h % J ~ 1 are 
the Yukawa coupling constants. 

The scale of neutrino masses (10 11 — 10 12 times smaller than the electron and quark 
masses) can be a result of the see-saw approach [3-5] . It leads after the integrating out of the 
corresponding heavy states (e.g. the right handed Majorano neutrinos with heavy masses of 
an order of M ~ 10 14 — 10 15 GeV) to the appearance of higher order effective Majorano mass 
operator (see [3], [5]): 

(ui(M u ) ij uj) = (tpl/M)^^), 

where h lj ~ 1 are the neutrino coupling constants and M v = {ip\/M)h is the neutrino mass 
matrix. 

2. Quarks mass matrices and the analytical form of the CKM matrix. 

The running of u- and d-quark masses calculated [7] in the first (dashed curves) and in the 
fourth (solid curves) order of QCD perturbation theory are shown [8] in Fig. 1. The quark 
running masses mi(/i) are related in MS renormalization scheme at the scale fx — Mi to their 
pole masses Mj by the well known relation: 

671 71 

where K=ll,2. The scale // = M t = 174.4 GeV is the most natural for the SUSY Standard 
Model. The running quark masses at this scale, (see in Fig. 1) have the following values (in 
GeV's): 

m u (M t ) = (0.21 ± 0.1) • 1(T 2 , m d (M t ) = (0.42 ±0.21) • 10~ 2 , m e = 0.511 • 10" 3 , 
m c (M t ) = 0.59 ± 0.07, m s (M t ) = 0.082 ± 0.041, = 105.66 • 10~ 3 , (7) 

m t {M t ) = 163 ± 4, m b (M t ) = 2.80 ± 0.40, m T = 1.777 ± 0.0003. 

These values differ strongly from the ones used in Ref . [2] which are determined at a very small 
scale ji — 1 GeV. However, the Stech's relation between them still holds. 
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Note that the matrices M d and M e in Eqs.(2),(3) have the block structure and their 



diagonalization is trivial. In fact, for any 2x2 matrix m 
contain in their left-upper part, one has: 

m diag = umvT = 

V0 ^ 

Here t 2 # = tan 2$ = 2p/(a — b) and s = sin-$ 



a p 
P b, 



which M d and M u 



u 



(J *), l*w = \(a + b±y/(a-by + 4f). 



(8) 
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-^(1 + (1 + 2 ) 2 ' anc ^ ^ ^ S the mixing angle. In Eqs.(2),(3) for the matrices M d and 



M e one has a = and the mixing angles i? e are small since t 2 # d = 2 ° , = 8cr < 1 and 
/3~o- 3 

is even smaller. Due to the block structure of M d and M e the unitary 
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matrices U d and U e also have the following block structure: 
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where s d = -±=(1 - (1 + (8a) 2 )~^ ~ 4<r(l - 24a 2 ) ~ 0.214, s e = - (1 + Qa 2 )-^)h ~ 



<t(1 



o- 



0.0705 up to the terms of the order of <r 4 . 



Let us note that Stech's matrices M d ,M e in Eqs.(2),(3) can be reconstructed using their 
diagonal form (i.e. the physical masses of d-quarks and electrons): 

Mf a9 = ( m -m s ) , Mf a9 = ( m -m, ) , 

by means of Eq.(4), which states: 

M d = UtMf a9 U d , M e = U + Mf a9 U e . 



(10) 



(11) 

The matrices U d = 6f 2 ( or U e = 6l 2 ) in Eq.(9) can be considered as rotating the 12 generations 
of d-quarks (or electrons). 

A bit more complicated (but much more instructive) is the diagonalization of the ma- 
trix M u . Similarly to M d and M e this matrix can be represented as : 



M u = U + M^U U , where U+ = 6+6+ (S) (6 U 12 )+ 



(12) 



Here the matrices: 

' c u -s u \ 

(6 u 12 )+ = s u c u , 6+ (6) 

{ 1/ 



C13 s 13 e- iS \ 
1 
v -s 13 e lS C13 / 
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rotates 12, 13 and 23 generations, respectively. The 13 rotation includes naturally the complex 
phase 5 which violates the CP-parity conservation of a theory. It can not be removed by a 
trivial phase transformation of the u-quark or t-quark fields. However, the value of s u turns 
out to be very small (s u ~ a 5 < 1) and one can really put 0" 2 ~ 1. The quark CKM matrix 
is: 

VSkm = U:U d = d+6+(5)6+, where 6+ = {6 u l2 ) + 6 d 12 ~ O d 12 . (14) 



Here in general (0" 2 ) + = 



Cl2 Sl2 



\ 



Sl2 C12 

V o o ij 

since s u ~ a 5 (see below) is neglegibly small. 
Thus for the given upper quark masses M dia9 



, with S12 = s d l2 c u -c d l2 s u = sin(tf d -^ 2 ) ~ sin$ d 



m c and the given quark mixing 
angles (from the quark CKM mixing matrix V^ KM) see below) the form of [/+ is fixed: 
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The same matrix £/+ can be obtained in the form of decomposition in powers of a (up to a 4 
terms) by a direct calculation: 
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Multiplying C/+ in this form by Ud one obtains with the same accuracy the following CKM 
quark matrix Vqkm = U^Ud'- 
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Comparing the CKM matrix (5) with Eq.(17) one finds: 

Sl2 = Sd ~ 4^(1 _ 24a 2 ) « 0.215, s 23 = -£=(1 - |a 2 ) « 0.0408, 



a 2 



implying that \R\ 



S13 = s 12 s 23 \R\ = ^ 2 (1 - ^ 2 /4) w 3.36 • 10- 3 
1 0.38 



(17) 



(18) 



S12S23 ^8(i _ 25a 2 ) 
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m 3 (19) 



Calculating M u by means of Eq.(12) one first obtains M' u = Of 3 M^ ia9 13 and then finds 
M u = 6i 3 M' u 6 23 in the form: 

f s n s 23 e- i5 s 13 e~ iS \ 

-Sl3 s 23 e S 23 + mi S23 

\ Si3e J<5 s 23 1 / 

up to the terms of an order of a 4 

Eq.(19) really reproduces the Stech's matrix M u in Eq.(2) only for negative sign of m 2 = 

o 2 

— |m 2 |; for positive m 2 = |m 2 | the ratio [i 22 — (M u ) 22 /m t is ^a 2 instead of — in Eq.(2). 
Considering higher order a 2 corrections to the matrix (19) one finds that /in = (M u )u/m t 
turns out to be 2a 4 for positive m u = \m u \ and is much smaller, — for negative lightest 
eigenvalue when m u = —\m u \. 

Therefore, literally the Stech's matrix M u is reproduced in the form (2) only for negative 

. . (-\m u \ 
mi and m 2 , e.g. at: M^ a9 = —\m c \ 

Thus, the u-quark mass matrix M u can be reproduced by the former of Eqs.(12) using 
the U u and C/+ matrices in the form of the latter of Eqs.(12) (or by Eq.(15) if (Oi 2 ) + — !)■ 
Actually this method will be very useful later for the restoration of the neutrino mass matrix. 

All these results are in a good agreement with the experimental data presented above. Also 
the diagonalisation (4) of M u and M d with the help of U u and U d matrices leads to the fol- 
lowing reasonable quark masses (in GeV's, as in Eq.(l)): 

m u = a\l - a 4 )m t (M t ) ss 2.00 • 10~ 3 
m c = tr(l - a 2 /2)m t (M t ) w 0.56 
m t = (1 - a 2 /2)m t {M t ) w m t {M t ) = 163 



m d = 1(^/1 + (8a) 2 - l)m b = 0.416 • 10" 3 



m s = + (8a) 2 + l)m b = 0.085 

m b = m b (M b ) = 2.83 

(20) 

These results correspond fairly well, as Stech has remarked [2], to the experimental data (7). 

The minus sign of some eigenvalues in Eq.(10) and in M^ mg can be easily removed by 
redefinition of the corresponding quark field: — * ^qu- 



3. Neutrino mass matrix, the leptonic CKM matrix and CP-phase 

To proceed further let us introduce the neutrino mass matrix M v with a structure similar 
to that of M u . The Super Kamiokande data [9] suggest a large v il y T neutrino mixing i.e. large 
t 23 = tg 2^3 > 1, or sin 2 2^3 = (1 +^ 23 2 ) _1 - 1 (i-e. s 23 ~ c 23 ^ ^ ^ 0.71) and not too 
small value of Am| 2 = rn 2 — m\ ~ (0.59 ± 0.20) 2 10" 2 eV. As m 3 3> m 2 (see below) one has: 

m 3 ~ y/Ami 2 ^ (0.059 ± 0.020)eV (21) 
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Simultaneously the atmosphere and solar neutrino data [10 - 12] (see the discussion in 
Refs. [13,14]) show a large suppression of u e u^ and also u e u T oscillations which have not yet 
been observed. This can be a result of small mixing angles si 2 < S13 <C S23 ~ 1/a/2 (see Ref. 
[10] and also [12 - 14]): 



s 12 ~ 0.035 ±0.020, si 3 ~ 0.25 ±0.10, s 23 = 0.7 ~ 1/V2 (22) 

and some (obviously not too large) hierarchy of three neutrino masses m-y <C m 2 <C m 3 of a 
type considered above for the quarks and elections: 

m 3 : m 2 : mi = 1 : al : ct*. (23) 

Here a"l is an unknown parameter which we can choose to be equal to a = 0.058 to avoid 
the introduction of additional new parameter a v = V0.058 ~ 0.24. This gives m 2 — 0.34 • 
10~ 2 eV and together with Eq.(22) seems to be in approximate agreement with the atmospheric 
neutrino data [12]. 

Essentially the other possibility has been suggested in a recent paper [15], where the 
neutrinos vi, v 2 and v 3 have been considered with almost equal masses m 3 ~ m 2 — mi ~ 3eV 
but a large hierarchy has been introduced in their mixing angles 

sis - < sia ~ s 23 ~ l/y/2 (24) 

This situation is not considered below as it seems more natural that neutrino have small 
mass hierarchy of the type (1) like the u-quarks, but with much smaller power of hierarchy: 
<j v ~ 0.24 3> <j. It is very suitable to describe the situation in neutrino physics since the 
condition (23) m 2 < m 3 (and mi < m 2 ) together with Eqs.(22) suppress strongly the non 
observed (at least for a moment) z/ e z/ M oscillations from electron neutrino sources on the Earth. 

Therefore, let us take approximately the central values of the data given in Eqs.(21)-(23) 
as a basis of our approach (see e.g. Ref. [5]) putting also a v ~ V0.058 in Eq.(23) and taking 
the following values for the neutrino masses: 

(mi \ /0.019 \ 

M dtag = ma = QM lO-^eV (25) 

V mj V 5.90/ 

where m 2 ~ 0" 2 m 3 , mi ~ <r 2 m 2 and the central value of m 3 has been taken from Eq.(21). 

The main point of this paper is that approach of the type of Eqs.(12)-(16) together with 
the electron matrix M e (determined by Eq.(3)) allows one to construct in general analytical 
form the neutrino mass matrix M v and also the leptonic CKM matrix. The matrix M v 
so obtained will have all desired properties reproducing naturally Eqs.(21) - (23) and will 
naturally include the CP-phase. 
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To this end let us remind that M e in Eq.(3) has the form: 

M e = U+Mf a9 U e , 



(26) 



where U e = O 



12 



/ c e s e \ 

S e C e 



has been determined in Eq.(9) with s e — 0.0705. Here 



M diag 



rrie 



m. 



rn-r 



V 001/ 

/0.5175 

103.6 
V 1777/ 



Y 



MeV represents exactly the three electron masses 0. 
The general form of the matrix M v can be written quite similar to Eq.(12): 

M v = U+M? a °U u with U+ = 6 l 2 td l +{5')dli , (27) 



where OX. 



12 



is defined similarly to 0\ 2 with the substitution s v = sini?^ 



/ c u s u \ 

-s v c v 
V 1/ 

for s e in 0\ 2 defined above (the value of i?5l 2 wm be determined later on). The matrices: 



( c 13 s 13 e~ l5 ' \ 
1 







O 



23 



Cl3 / 



( 1 M 

C23 S23 
V -s 23 c 23 / 



f28) 



represent the rotation of 13 and 23 generations of M^ m9 , respectively, and 5' is the leptonic 
CP-phase which has appeared here naturally (clearly Sij, Cij from here and later on means the 
sines and cosines of the leptonic and not quark mixing angles). 

To calculate M v in Eq.(27) in explicit form we note that according to Eq.(8) one has 



w 12 



'mi 



m 2 



O 



m 3 y 



12 



/ a p \ 

p b 
V m 3 / 



(29) 



with a = s 2 v m 2 + c 2 v rti\ ~ mi, b = c 2 v m 2 + s 2 v m\ ~ m 2 , p = —c u s v {m 2 — mi) and \p\ <^ m 2 is 
small since s u in Eq.(22) is very small (actually s v ~ si 2 — s e — 0.036 see below). Then it is 



1 Let us mention that the matrix M e can be slightly modified: M e 



ML 



( Xa 4 yfa 3 
fa 3 -/3a 2 



a 



V a / 

where /3 = 1 + Ai with very small Ai = 0.0202308 and A = 0.0106891 adjusted to reproduce exactly the well 

/0.510999 \ 

I 



known masses of all three electrons: M e dla9 — U e M e U+ 



I MeV. This modification 



105.6584 
V 1777.0/ 

will change negligibly the leptonic CKM matrix, leading to s e ~ 0.0689 and is not important at all as the 
neutrino basic parameters are determined very roughly in Eqs.(21)-(23) (see also the Conclusion). 
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easy to calculate the matrix: 



m> = ous') 



I a p \ 

p b 
V m 3 J 



6is05') 



and further to find the following result for the neutrino mass matrix M v = Oi^M'O 



i . 

23- 



M u = 



ci 3 si 3 S23(^3 - a)e tS ' 
+pci 3 c 23 
c 2 3Ci 3 si3(^3 - a)e l<5 ' 
V -PS23C13 



Ci 3 Si 3 S23(^3 - a)e- lS ' 

+PC\ 3 C 23 

S2 3 (^3C? 3 + a 2 s? 3 ) 

+C2 3 & + A 22 

c 2 3S23("^3C? 3 -b + as\ 3 ) 

+ A32 



C23Ci3Si 3 (^3 - a)e~ l5 ' \ 

-ps 23 c 13 

C 2 3S23("^3C2 3 ~ b + asj 3 ) 



-A 



23 



c 23( m 3C? 3 + as 



13 J 



(30) 



+sl 3 6 + A 33 J 



The values of A 2 2 = — psi 3 sin 2$2 3 cos 5', A 23 = — ps\ 3 (e~ l& ' — 2s2 3 cos5') = A 32 and of 
A 33 = — A 22 are small and can be neglected. Also since a <C 5 <C m 3 , s u ~ s i2 (see below) 
and s\ 2 <C Si 3 <C 1 are very small, one can disregard in the matrix M v all the terms containing 
a, b and put c i2 — c i3 ~ 1. The value of |p| ~ Si2"^2 is also very small \p\ <C Si 3 m 3 , but the 
terms containing it in the matrix (30) can not be omitted as that would violate the normal 
complex structure of the matrix M v and of CKM matrix considered below. Therefore omitting 
small terms one obtains the matrix M v in the following simple form: 



My = 



s\ 3 m 3 + ac\ 3 



ci 3 (s 13 s 23 m 3 e i5 ' + pc 23 ) Ci 3 (si 3 c 23 m 3 e iS ' - ps 23 ) \ 



ci 3 {s 13 s 23 m 3 e- +pc 23 ) 
V c 13 (s 13 c 23 m 3 e tS ' - ps 23 ) 



s 2 23 m 3 + c\ 3 m 2 
S23 C 23( m 3c| 3 - m 2 ) 



s 23 c 23( m 3 c 13 



m 2 ) 



4 3 m 3 + s\ 3 m 2 



(31) 



Here, s 2 3 ^ c 23 ~ l/V~2 and not too small values of s± 3 are determined in Eqs.(22). 
In conclusion of this section let us construct the leptonic CKM matrix: 



v 1 ckm = u^Ue = d l +d[t(6')d[ 2 , 



(32) 



where 



w 12 — w 12 w 12 — 



V 



C12 S12 N 

-S12 C12 
001; 

and S12 = sin(tff 2 - ^12) = sin^ 12 ~ 0.035 ± 0.020 as it is determined by Eq.(22). This gives 
for the neutrino v x v 2 mixing angle: i?5l 2 = ®\2~ $12 — 0.036±0.020 (in radians, or (2.1±1.1)°). 
Multiplying the matrices O l £ in Eq (32) one obtains Vq KM just in the form of Eq.(5) with 
S23 - c 23 - and s 12 , s 13 given by Eq.(22). 
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Let us emphasize that the approach of this section, in particular the basic parameters 
in (21), (22), and also the relation (23), are approximate and have a model nature: also 
the hierarchy parameter o v is the main one and its value a v = 0.24 is determined only 
approximately by modern experimental data. 

Numerically the neutrino mass matrix (31) is represented in our approach by: 



M, 



/ 0.130 0.335e~ i<5 ' + 0.27 • 10~ 2 0.341 e - i<5 ' + 0.26 • 10~ 2 \ 

0.335e i(5 ' + 0.27 • 10~ 2 0.965 0.868 

v 0.3Ale i5 ' - 0.26 • 10~ 2 0.868 1 



and the CKM leptonic matrix (32) is 



V 1 

V CKM 



( 0.968 0.0339 

-0.175e i(5 ' - 0.025 0.714 - 0.061e i<5 ' 
v -0.178e i<5 ' + 0.025 -0.699 - 0.0863e i<5 ' 



0.25e~ i<5 ' \ 
0.678 
0.692 j 



m 3 /1.975 
(33) 

(34) 



Unfortunately instead of predicting the neutrino mixing angles Si 2 , s± 3 , s 23 we have used 
their values (22) which are badly defined by experiment. The hierarchy parameter a v remains, 
as has been mentioned above, practically free and we have put it to ^/o : ~ 0.24 just by hand. 
Both the matrices 0\ 2 and 0{ 2 represent a simple Abelian rotation: 0\ 2 by the angle $ v l2) 0\ 2 
by the angle (— i?^) an d ^12 by the angle d\ 2 . Therefore the product 0\ 2 d{ 2 = 0\ 2 0\t leads 
to a rotation by the angle d l 12 = d\ 2 — i?i 2 , where $\ 2 ~ 2§\ 2 (the value of s e has been given 
above just after Eq.(9) and s 12 in Eq.(18)) and therefore -Q\ 2 ~ -&[ 2 , or s 12 ^ s u . 

Similarly to the case of the quark CKM matrix, the most natural value of the leptonic 
CP-phase leading to the largest possible CP violation can be 5' = tt/2 or rj' = e l5 ' = i. This 
CP-phase 5' can manifest itself in the Pontecorvo neutrino oscillation experiments. It is very 
difficult to observe it now. Below we discuss shortly the possibility of these observations. 
The exact expressions for probabilities of neutrino oscillations are given in Appendix, since 
they are very cumbersome. Some of them have been obtained earlier in a number of papers 
[16 - 20]. 



4. The leptonic CP-phase in neutrino oscillations experiments 

Many papers were devoted to the studies, pioneered by Bruno Pontecorvo [16], of two and 
of three [17 - 20], [5] neutrino oscillations. We consider them below shortly in order mainly to 
specify the role of the leptonic CP-phase [17] in these oscillations. 

Let us express v e , and v T fields entering the weak interaction Lagrangian in terms of 
neutrino states vi, z/ 2 , z/ 3 with definite masses mi,m 2 ,m 3 using the leptonic CKM matrix (5) 
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(or (Al) from Appendix) as follows: 



u e (ct) = c 13 z/i(0)e- i£l * + si 2 ci 3 ^(0)e- i£2 * + s 13 i/ 3 (0)e-* £3 '- i<5 ' 

u^ct) = -{s 12 c 23 + s 13 s 23 e iS ')v 1 (0)e-^ t + c 23 z/ 2 (0) e -^* + c 13 s 23 i/ 3 (0) e - i£ 3* (35) 
Md) = (s 12 s 23 - s 13 c 23 e l5 ')u 1 (0)e^ t - s 23 z/ 2 (0)e^ 2 ' + c 13 c 23 z/ 3 (0)e- l£3 ' 

up to the terms of the second order in small quantities Si 2 ,Si 3 <C 1. (see Appendix for the 
exact Vqkm matrix) Here t = L/c is the time when neutrinos are observed at a distance 
L = ct from their source; the probabilities of neutrino observation at this distance from the 
source is P a b(t) = |(z/ a (t)z/ fc (0))| 2 . Multiplying u a (t) from Eq.(35) by P&(0) and taking into 
account the orthogonality of u a (Q) states (v a (0)vb(0)) = 5 a b, it is easy to find: 

P{y e v e ) = |c 2 3 + s 2 2 c 2 3 e^ + s 2 3 e^ 1 | 2 

P(W) = ||si3C 2 3 + si3S 2 3e l5 '| 2 + c2 3 e^+c2 3S 2 3 e^3i|2 ( 36 ) 

P{v T V T ) = ||Sl 2 S 2 3-S 13 C 2 3e i5 '| 2 + s2 3e ^l +C 2 3C 2 3e ^3l|2 



P{y e vp) = \c 13 (c 23 s 12 + s 23 s 13 e iS ') - s X2 c X3 c 23 e^ - s 13 c 13 s 23 e^ : ' + ^\ 2 
P{v e v T ) = \c 13 (s 23 s 12 - c 23 s 13 e i5 ') - s 12 c 13 s 23 e-^ 21 + s 13 c 13 c 23 e t(5 ' + ^\ 2 
Piy^r) = I (s 12 s 23 - s 13 c 12 c 23 e~ t5 ')(s 12 c 23 + s 13 c 12 s 23 e tS ') + c 23 s 23 e lLp21 

-s 2 3C 2 3C 2 3 e^ 31 | 2 



(37) 



where 



and 



= <« - OW' = ^ = (38) 

Am 2 - = m 2 - m 2 , i,j = 1, 2, 3 (as Si ~ + —-^ at cp > m^. 

Eqs.(36)-(38) determine the neutrino oscillation probabilities in vacuum ignoring the very 
important in some cases MSW effect of the medium influence (see in ref.[10]). This effect 
has been well studied and can be included separately (e.g. for v e v e , u e u^, u e u T in the solar 
neutrino case). The expressions for P(v^v e ), P(v T is e ), P(y T v^) coincide with Eqs.(37) with 
the substitution 5' — > —5' 

The simple algebra allows one to reduce Eqs.(36),(37) to the partially known [17 - 20] 
expressions containing the leptonic CP-phase 8' and given below in Appendix. 

Let us rewrite the Eqs.(36),(37) in much more simple forms calculating numerically the 
coefficient in front of cos ipij = 1 — 2 sin 2 ((^j.,/2) using for example the central values of leptonic 
mixing angles -d\ 2 , i?' 13 , # 23 given in Eqs.(22), i.e. their sines si 2 ,si 3 ,s 23 , respectively. These 
equations in their numerical form show clearly the influence of the leptonic CP-phase 5' on 
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the neutrino oscillations patterns: 

' 1 - P{v e v e ) ~ c\ 2 sin 2 (2< 3 ) sin 2 (^ 3 i/2) + cf 3 sin 2 (2^ 12 ) sin 2 (y9 21 /2) 

+s 2 2 sin 2 (2^ 13 )sin 2 (v932/2) ~ 
0.23 sin 2 (v?3i/2) + 4.3 • 10~ 3 sin 2 (y9 2 i/2) + 0.29 • 10~ 3 sin 2 (^ 32 /2) 
< 1 - P[y^) ~ (0.94 + 0.017 cos 5') sin 2 (y9 32 /2) + (0.062 - 0.017 cos 5') sin 2 (v? 3 i/2) (39) 

(0.063 + 0.016 cos 5' - 0.31 • 10~ 3 cos 2 5') sin 2 (y9 21 /2) 
\-P{y T v T ) ~ (0.94 + 0.017 cos 5') sin 2 (v? 3 2/2) + (0.062 - 0.017 cos 5') sin 2 (y9 31 /2) 

(0.064 - 0.016 cos 5' - 0.31 • 10~ 3 cos 2 5') sin 2 (y9 21 /2) 

Eqs.(37) determine the probabilities for neutrino of a given sort z/j to change its type (i.e. 
to transform into the other sort Vj) at a distance L from the source of neutrinos. For v e 
this probability does not depend on the CP-phase 5' at all as can be seen already from 
Eqs.(36),(37), however, for transformations of and v T this dependence is not too small. 

The probabilities for different neutrino transitions P{yiu^) ^ P{yjV^) depend also on the 
value of 5': 

' P(z/ e z/ M ) ~ 0.059{1 - 0.14 cos 5' - 0.98[cos(y? 3 i) + 0.019 cos(y? 2 i) + 0.0012 cos(y? 32 )] 

+0.140[cos(v? 32 + 5') - cos(y? 31 + 5') - cos(y? 2 i - 5')]} 
P{y e v T ) ~ 0.061{1 - 0.13 cos 5' - 0.98[cos(y? 31 ) + 0.018 cos(y? 21 ) + 0.0012 cos(v? 32 )] 

-0.135[cos(y? 31 + 5') + cos(y? 32 + 5') - cos(y9 21 -5')}} (40) 
Piy^r) ~ 0.47{1 - 0.99[cos(y? 32 ) - 0.066 cos(v? 2 i) + 0.061 cos(^ 3 i)] 

+0.035 sin 5' cos( ^ 31 ^ 32 ) + 0.035 sin 5' cos ^ 31 ^ 32 ) sin(y9 21 /2) 

-0.017 sin <S'sin(^ 2 i)} 

It is interesting to consider the time reversal effect which reveals in the difference between 
the probabilities P[y e Vy) and P(f M z/ e ) or P{v e v T ) and P(v T v e ) etc. [17]. Taking the difference 
between P{y e v^) — P(z/^z/ e ) and the same for v e v T , v^Vr one obtains: 

P(V(iVe) - P(v e Vn) = 0.0164(siny? 21 + sin^ 32 - siny? 31 ) sin 5' 
< P(v T v e ) — P{y e v T ) = — 0.0164(sin ip 2 i + sin y? 32 — sin </? 3 i) sin 5' (41) 

P(^) - P(u^r) = 0.0328(cos - cos (v?31 + V?32) ) sin<5' 

where 0.0328 = Ci 3 sin 2$' 12 sin 2$' 13 sin 2i? 23 for the numerical coefficients given above. So, 
the leptonic CP-phase can manifest itself in time reversal neutrino transitions v a Vb and v^Va 
experiments. However, these experiments need the neutrino beams with a fixed value of 
Lj1E v what is very difficult to organize because usually one deals with continuum spectra of 
produced neutrino. 

We emphasize once more that the numerical coefficients in Eqs.(39),(41) can be determined 
by future experimental data only, but the general form of them (obtained algebraically from 
exact Eqs.(A2)-(A6) of Appendix) remains always valid. 
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Let us average Eqs.(39),(40) over L/2E U considering the case of large L/2E U > (Am^)" 1 , 
i.e. of large angles = Am]jL/2E u . This nicely corresponds to the real experimental 
situation where one works with a continuum spectra of neutrino. This averaging gives: 

(1 - P(v e v e )) ~ 0.1193, 
(1 - P(*VV)) - 0.529(1 + 0.016 cos 8' - 0.289 • lO^cosV), (42) 
(1 - P{v T v T )) ~ 0.531(1 - 0.015 cos 8' - 0.288 • lO" 3 cos 2 8') 



and 



(P(z/ e z/ M )) ~ 0.0585(1 + 0.140 cos 5'), (P(^r)) ~ 0.0608(1 - 0.134cos5'), 
(P(Wt)) ^ 0.470(1 + 0.394 • 10~ 3 cos 8' - 0.163 • 10~ 3 cos 25') 



(43) 



where brackets mean the averaging over all ipij. 

Comparing Eqs.(39) and (40) (and also (A2)-(A6) in the Appendix) one finds: P{yiV^) + 
P(uiVj) + PiviV k ) = 1 for different v h Vj and v k (e.g. 1 - P{u e u e ) = P{y e v^) + P[y e v T )). 

Eqs.(42)-(43) show also that the leptonic CP-phase 8' can be observed experimentally, in 
principle, by measuring the average v e v^ or u e u T transition rates with 14% accuracy which is 
much larger than the effect (41) of time reversal. 

This is illustrated in Figs. 2(a,b) where the averaged probabilities of v^Vj transitions 
(P(viVj)) are calculated for the values of si2,si3,S23 determined in Eqs.(22). These figures 
show a large (about 30%) difference between {P(v^v e )) and {P(v T v e )) probabilities dependen- 
cies on the CP-phase 5'. For example the value of {P(v^v e )) ~ 0.067 at 5' = turns out to 

be larger than (P(u T u e )) ~ 0.057 by 17%, i.e. |p|^ £ jj - I- 17 , wnile at 5' = tt vice versa 

(P(v»Ve)) - 0.050 becomes smaller than (P{u T u e )) ~ 0.059 by 15% : 7W^4\ - 0.85. 

Unfortunately, the absolute values of these probabilities are small of about 1/20, but nev- 
ertheless they be really observed experimentally. Even weaker is the 8' dependence of the 
average value of the probability of v^v T transition rate (averaged over the oscillations con- 
nected with different L/E u values, or over (p^ - as in Figs. 2(a,b)), shown in Fig. 2c. This 
figure shows that (P(i/ j[1 i/ T )) changes by 0.05% only when 8' changes from to n. 



5. Conclusion 

The following three simple problems were discussed and solved in this paper: 

a) The u- and d-quarks and electrons mass matrices were obtained in a simple hierarchical 
form, including quark's CP-phase 8 and correcting the matrices suggested by B.Stech, 

b) The recent data on neutrino masses and mixing angles were discussed shortly and used for 
construction of the neutrino and leptonic mass matrices and CKM matrix, both including the 
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leptonic CP-phase 5', 

c) The CKM matrix obtained was used to investigate the three neutrino oscillations in the vac- 
uum. The method of determination of the leptonic CP-phase from their observation (averaged 
over the oscillations) was presented. 

Note, that there is a vast ambiguity in the determination of quark's (or leptonic) mass 
matrices. E.g. the pair of matrices M' u = U M d U^, with any unitary matrix U , leads exactly 
to the same mass eigenvalues and to the same mixing angles as M u , 

Also the following correction has to be added to the central part of the paper. It was noted 
there (in Section 2, after Eq. (21)) that Stech's matrices (2), (3) have negative value of some 
masses and that they positivity can be restored by a simple 75 transformation. However this 
will violate the symmetry of different generations fields and also will change the form of mass 
matrices. Better is to avoid this shortcoming and construct the particles mass matrices with 



only positive eigenvalues Mf a9 



>m 



<ii 



m 



directly from Eqs. (11) [or (29) and (12)], 



e.g. for d-quark and electrons one finds: 



K 



I a p \ 

p b 
V m b ) 



where < 



a 
b 



c\ 2 m 2 ■ 



s\ 2 m\ 



c 2 12 m 2 



M' d 



15a 4 



7 rr 3 



7 rr 3 



V 







V 










a 



P = -Ci2Si 2 [m 2 - mi J ~ 
/ 

and similarly: 



1 2 

2<7 TUb 



or: 



J 



V 



3A' l( 7 4 

I- 3 




(1 + P')a 2 









a 



J 



with small A' = 0.020, (3' = 0.016. Both M d , M e differs slightly from Stech's forms given in 
Section 1. The u-quark mass matrix is determined by Eqs. (12), (19) and also slightly deviates 
from the Stech's form given by Eq. (2): 



/ 2a 4 1 



\ 



a 2 i] 



if 



T2 



a 2 rj + \ 
1 / 



--a 



u \ + 



m t , where lit = O^Ot 3 (5)(O u 12 ) 



The diagonalization of all these matrices M^ m9 = U u M' a U+ can be done by the same unitary 
matrices Ud, U u , U e , U v with the same mixing angles (used at they construction) as were 
used above for the Stech's matrix case. 
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6. Appendix 

We give below the probabilities rates PiviVj) in general algebraically form using the well 
known exact leptonic CKM matrix (given above for quarks in Eq.(5)): 



i 

CKM 



( C12C13 
— s 12Cl3 

V S12S23 - ci 2 c 2 3Si 3 e 



Cl2-S 2 3' s 13 e 
iS' 



iS' 



S12C13 
C12C23 — •Si2S 2 3Si3e 



iS' 



-C12S23 - Sl2C 2 3Sl 3 e' 



S23 c 13 

C23C13 J 



(Al) 



At 5' ^ it leads to: 



1 " PWVy) 



{4 sin 2 (2tf 21 ) + 44 sin 2 (2tf 23 ) + 4»4» sin 2 (2< 2 ) + 44 sin 2 (2tf 



23; 



+ cos5 , sin(4^ 12 )sin 2 (2^ 23 )(si34 - 44) - cos 2 5's 2 l3 sin 2 (2tf 23 ) sin 2 (2^ 12 )} sm 2 (</? 21 /2) 



+{44sin 2 (2tf 23 ; 
+{c 2 2 c 2 3sin 2 (2^ 23 ; 



+ 44 sin 2 (2< 3 ) + cos5'sl 3 ci 3 sin(2^ 12 ) sin(2tf 23 ) sin(2tf' 13 )} sin 2 (y9 31 /2) 
+ 44 sin 2 (2^ 13 ) - cos5's| 3 Ci 3 sin(2^ 12 ) sin(2tf 23 ) sin(2^ 13 )} sin 2 (y9 31 /2) 
' " ' (A2) 



1 - P(v T v T ) = {4 sin 2 (2tf 21 ) + 44 sin 2 (2tf 23 ) + 44 sin(2< 2 ) + 44 sin 2 (2tf 23 ) 

+ cos 5' sin(4^ 12 ) sin 2 (2^ 23 )(44 - s 13 4) - cos 2 5's 2 13 sin 2 (2tf 23 ) sin 2 (2^ 12 )} sin 2 (y9 21 /2) 
+{44 sin 2 (2tf 23 ) + c 2 2 4 sin 2 (2< 3 ) - cos 5'c 2 3 c 13 sin(2^ 12 ) sin(2tf 23 ) sm(2<4)} sin 2 (y9 31 /2) 
+{ 44 S ^ n (2^ 3 ) + 44 sm 

(2^ 13 ) + cos<5'4c 13 sin(2< 2 ) sin(2tf 23 ) sin(2^ 13 )} sin% 32 /2) 

(A3) 

for 1 — P(v e is e ) see Eqs.(39). As before here (pij = Am 2 j L/2E l/ and: 

P(y*Vp) = |{sin 2 (2^ 13 )(4 + 44 + 44 + \c l3 sin(2tf 23 ) sin(4^ 12 ) cos 5') 

-2c 2 3 sm 2 (2^ 12 )(4 - 44) cos(v?2i) - 24 sin 2 (2^ 13 )(c 2 2 cos(v? 3 i) + 4 cos(y9 32 )) 
+C13 sin(2i?' 12 ) sin(2i? z 13 ) sin(2i? 23 ) (4 cos(<5' + tp 2 i) — c\ 2 cos(5' — v?2i)) 
+C13 sin(2^ 12 ) sin(2tf' 13 ) sin(2tf 23 ) (cos(5' + y? 32 ) - cos(5' - </? 3 i)) 

+244sin 2 (2< 2 )} 

(A4) 

P(y e v T ) = |{sin 2 (2^' 13 )(4 + c\ 2 c\ 3 + 44 - \c 13 sin(2tf 23 ) sin(4^' 12 ) cos 5') 

+24 sin 2 (2^ 12 )(4 - 44) cos(v?2i) - 24 sin 2 (2i? i 13 )(c? 2 cos(y? 3 i) + 4cos(</? 32 )) 
+C13 sin(2i^ 12 ) sin(2i? J 13 ) sin(2i? 23 )(4 cos(<5' — </? 2 i) — s\ 2 cos(<5' + V21)) 
+C13 sin(2tf' 12 ) sin(2tfy sin(2tf 23 )(cos(<5' + </? 3 i) - cos(5' + ^32)) 

+244sin 2 (2^ 12 )} 

(A5) 
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\ {2sj 3 s in 2 (2^ 12 ) cos 2 (2tf 23 ) + {c\ 3 + c\ 2 + s\ 2 + (cf 2 + s\ 2 )s\ 3 ) sin 2 (2tf 



1 ') 

23/ 



-[2s 



2 (A 



4 3 )sin 2 (2^ 12 ) + [2 S 2 3 (4 + 4) 



■[2c 



13^12 



(1 + sj 3 ) sin 2 (2^ 12 )] sin 2 (2tf 23 )] cos(y? 21 ) 
c\ 2 s\ 3 ) sin 2 (2$ 23 ) - ici 3 sin(2^ 12 ) sin(2^ 13 ) sin(4$ 23 ) cos 5'} cos(y? 3 i) 
+ [2c 2 3 sin 2 (2tf 23 ) (s\ 2 s\ 2 - c\ 2 ) - |c 13 sin(2^ 12 ) sin(2^ 13 ) sin(4tf 23 ) cos 5'] cos(y9 32 ) 
+2ci 3 sin(2# 12 ) sin(2^ 13 ) sin(2tf 23 ) sin 5' sin(^ 2 i/2) cos( ^ 1 + ^ 32 ' 



+si 3 sin(4tf' 12 ) sin(4$ 23 ) cos5'[l + s 2 13 ] sin 2 (y9 2 i/2) 
— C13 sin^^) sin^-i^g sin(2i9 23 ) sin <5' sm(</? 2 i) 
-2s 2 3 sin(2^ 12 ) sin(2$ 23 ) cos(25') sin 2 (y9 21 /2)} 

(^46) 

For the average values of these probabilities over all ipij (i.e. over E u , or over L at large 
ifiij) one obtains the values the (42), (43) slightly dependent on the leptonic CP-phase 5': 



f (1 




= A ee ~ 0.1193 


(1 




= i w + B m cos <5' + C w cos 2 5' 


< < x 


- P{v T v T )) 


= A TT + B TT cos 5' + C rr cos 2 5' 






= A e)l + B efl cos 5' 






= A eT + B eT cos <5' 






= A^ T + 5 MT cos 5' + C MT cos(2<5') 



(A7) 



where the coefficients are: 



A ee — 
A^ = 



'■efi 



.2 

'13 1 
„2 



sin' 



2^13 



+ (C? 2 + S 



'12J S 13 



! (2^ 13 )] 
sin 2 (2$f 



23/ 
4 



A TT — 



+ ( S f 3 sin 2 (2< 2 ) + sin 2 (2^ 13 )) S23 

+4sin 2 (2^ 12 )] 

|[(c 2 3 + (4 + ^) S 2 3 )sin 2 (2^ 23 ) 

+ ( S f 3 sin 2 (2^ 12 ) + sin 2 (2^ 13 ))4 

+4 3 sin 2 (2^ 12 )] 
I[(l + C 4 2 + S 4 2)s 2 3sin 2 (2 ^ 3)] 

+2c 2 3 c 2 3 sin 2 (2^ 12 )] 
A er = i[(l + 4 + 4)c 2 3 sin 2 (2^ 13 )] 

+2c 2 3S 2 3 sin 2 (2^ 12 )] 
V = i[2s 2 3sin 2 (2^ 12 )cos 2 (2^ 23 ) 

+ sin 2 (2# 23 ){(4 + 4 2 ) S 2 3 

+ C 13 + C 12 + S 12}] 



B 



B 



C, 



h 13 sin(2tf 23 )( S 2 3 s 2 3 - c 2 3 ) sin(4tf< 



12/ 



C m = -is 1 3sin(2^ 23 ) Sl3 sin(2^ 23 )sin 2 (2< 2 ) 



B TT = 
C TT = 



i«i3 sin(2^ 23 )(s 2 3 - c 2 3 s 2 3 ) sin(4^ 12 ) 
■isia sin(2^ 23 )s 13 sin(2^ 23 ) sin 2 (2^ 12 ) 



|c 1 3sin(2^ 13 )sin(2^ 23 )sin(4^ 12 ) 

-|ci 3 sin(2#' 13 ) sin(2tf 23 ) sin(4^ 12 ) 

|(l + S 2 3 ) Sl 3sin(4^ 12 )sin(4^ 23 ) 
-| S 2 3 sin 2 (2^ 12 )sin 2 (2^ 23 ) 



(A8) 



The numerical value of these coefficient are given above in the text in (42), (43) for neutrino 
mixing angles Si 2 , S13, S23 values determined in Eqs.(22). 
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Figures captions. 



Fig.l. a) The run of the upper quark masses m t ,m c ,m u calculated in Refs.[8] in the 
first order of QCD perturbation theory (the dashed lines) and in the fourth order 
of it (solid lines). The vertical line corresponds to the scale /i = M t ~ 174.4GeV 
used in the paper. The run of electron's masses m r ,m M ,m e is disregarded (while 
it can be easy taken into account and is not essential), 
b) The same for the masses m^m^m^ of the lower quarks. 

Fig. 2. a) Dependence of the average (over oscillations) value (P(z/ e z/ M )) of z/ e z/ M transi- 
tion probability on the CP-phase 5'; at 5' = ir it has about 13% minimum. 

b) The same for the average value (P{v e v T )) of v e v T transition probability. As is 
seen instead of the minimum at 5' = ir, as was for the z/ e z/ M transition case, it has 
here the 14% maximum. 

c) The same dependence on the CP-phase as was shown in the cases a), b) but 
for the average v^v T probability (P(z/ M z/ r )). Its dependence on CP-phase 8' here is 
much more flat - about in hundred times smaller then in the z/ e z/ M and v e v T cases. 



20 



Fig.l(a,b) 
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Fig.2(a,b,c) 
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